ResearchGate 


See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/356915462 


Critical Angle and Fundamental Frequency of Symmetric Airfoils at Low 
Reynolds Numbers 


Preprint - May 2022 


CITATIONS READS 
0 30 
1 author: 

—  DilekFunda Kurtulus 


Neo 
x» Middle East Technical University 


158 PUBLICATIONS 719 CITATIONS 


SEE PROFILE 


Some of the authors of this publication are also working on these related projects: 


Project Bioinspired flapping wings View project 


Project Verification for Experimental Unsteady Aerodynamic Loads of a Flexible Flapping Wing Actuated with PZT Material View project 


All content following this page was uploaded by Dilek Funda Kurtulus on 10 December 2021. 


The user has requested enhancement of the downloaded file. 


Journal of Applied Fluid Mechanics, Vol. 15, No. 3, pp. x-x, 2022. 
Available online at www.jafmonline.net, ISSN 1735-3645, EISSN 1735-3645. 


Critical Angle and Fundamental Frequency of 
Symmetric Airfoils at Low Reynolds Numbers 


D. F. Kurtulus" 


! Prof. Dr., METU Aerospace Engineering Department, 06800, Ankara, Turkey 
+ Corresponding Author Email: kurtulus @ metu.edu.tr 


ABSTRACT 


The unsteady vortex evolution of NACA 0012 airfoil is numerically investigated at low angles of attack 
ranging from 0° to 10° where the separation is performed from the trailing edge. The Reynolds number ranges 
between 1000 and 4000. The laminar separation bubble at the trailing edge is observed and the main flow 
features are presented. It is found that, the increase of the angle of attack and Reynolds number result higher 
lift to drag ratio by an extensive decrease of the drag coefficient below 8? angle of attack. The transition from 
the steady condition to periodic force evolution has been revealed with a detailed flow field analysis at 
different angles of attack and Reynolds numbers. The critical angle is defined as the angle of attack where the 
onset of oscillations starts with a dominant fundamental frequency of oscillation. The angle of attack where 
the laminar separation bubble (LSB) is observed is also revealed in the current study. Both the LSB formation 
angle and the critical angle of attack is found to decrease with the increase of the Reynolds number from 1000 


to 4000. 


Keywords: Micro air vehicles; vortex shedding; unsteady aerodynamics 


NOMENCLATURE 
C chord [m] po free-stream static pressure [Pa] 
C; skin friction coefficient, tw/q q dynamic pressure (1/2pUo?) [Pa] 
Cp pressure coefficient, (p-po)/q Re Reynolds number 
Ci instantaneous lift coefficient, L'/qc st Strouhal number (fc/Uo) 
C, mean lift coefficient Uo — freestream velocity [m/s] 
Ca instantaneous drag coefficient, D'/qc V total velocity vector [m/s] 


- mean drag coefficient 
D' drag per unit span [N/m] 
dt time step [s] 


for fundamental frequency at acr [Hz] 
L' lift per unit span [N/m] 
p static pressure [Pa] 


1. INTRODUCTION 


The advance of technologies increases the 
interest to micro and even nano air vehicles which 
are imitating the flights of birds and insects. 
Reynolds number regimes for birds are in the range 
of 10* to 10? for birds and lower than 10^ for insects 
(Traub et al. 2019). Reynolds numbers of different 
insects which have also hovering capability are 
order of 10-10*. For example, fruit fly has a 
Reynolds number as low as ~150 in hover, 
similarly honeybee has Re~1200, dragonfly has 
Re~2000, hawkmoth has Re~6000 (Han et al., 
2014). 

The low Reynolds number flows also take the 
interest due to the unsteady aerodynamics around 


a angle of attack [°] 
Qcr critical angle of attack [°] 
aLss LSB formation angle of attack [°] 


At time increment [s] 
p fluid density [kg/m?] 
v kinematic viscosity [m?/s] 


the wind turbine blades. DNS simulations with 
spectral/hp element method are used by Nakhchi et 
al. (2021) to investigate the transition to turbulence 
and the laminar separation bubble around 
NACA4412 airfoil as a wind turbine blade airfoil in 
the Reynolds number range of 2.5x10* - 1.5x10° 
and the angle of attack range of 0°-16°. 

The airfoil performance is highly affected by 
the unsteadiness of the flow behind the airfoil (Lee 
et al. 1998). At low angles of attack, the separation 
is observed to start from the trailing edge of the 
airfoil at low Reynolds numbers which is also 
denoted by Huang et al. (2001) A laminar 
separation bubble (LSB) is observed at the suction 
surface of the airfoil (Jones, 1933). The 
reattachment point of the LSB is found to be 
unsteady which is also denoted to be one of the 
reasons of the shedding of the vortical structures 
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from LSB (Watmuff, 1999; Traub, 2008, Eljack et 
al. 2020, 2021). The investigation of the vortical 
patterns also reveals the high amplitudes 
oscillations in aerodynamic forces at these low 
Reynolds numbers. 

Ohtake et al. (2007) measured the aerodynamic 
forces and pitching moment with respect to quarter 
chord point of NACA 0012 airfoil in a 0.3 m x 0.3 
m cross section wind tunnel at Reynolds numbers in 
the range of 1x 10^ to 1x 10°. Counsil et al. (2013) 
have done a comparative study of two airfoils 
(NACA 0012 and SD 7003) at low Reynolds 
numbers ranging from 4.8 x 10* to 2.5 x 10° and at 
angles of attack ranging from O° to 8° using 
unsteady RANS shear-stress transport y-Ree model 
with different turbulence intensities and the ANSYS 
CFX solver. They found that by the angle of attack 
and Reynolds number increase, the laminar- 
separation bubbles shrank and moves toward the 
leading edge. Council et al. (2013) also compared 
their results and found similar results with the 
experimental data of Ohtake et al. (2007). Cleaver 
et al. (2012) perform experiments on a NACA 0012 
airfoil of dimensions 0.1 m chord and 0.3 m span in 
a free-surface closed-loop water tunnel with a 
turbulence intensity less than 0.5%. The lift force 
measurements for the stationary 2D NACA 0012 
airfoil are presented for Re = 10,000, 20,000, and 
30,000. They denoted that the lift curves for Re > 
20,000 are typical of leading-edge stall, and those 
for Re « 20,000 are typical of laminar trailing-edge 
stall (Cleaver et al., 2012). For Re « 20,000, as a 
result of the trailing-edge stall, the peak of the lift 
coefficient is more rounded and diminishes slightly. 
This phenomena is also observed by Huang et al. 
(1995) who have experimentally investigated the 
flow patterns and characteristics of vortex shedding 
and shear-layer instability of a NACA 0012 
cantilever wing using smoke wire and surface oil 
flow techniques. They have also measured the 
frequencies of the shed vortices in the wake region 
and the instability waves developed on the 
separated shear layers by using two one-component 
hot-wire anemometers (Huang et al.,1995). Kurtulus 
(2015, 2016) using laminar unsteady RANS 
simulations showed wake patters at Re=1000 for 
steady uniform flows around NACA 0012 and 
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Fig. 1 


NACA 0002 airfoils at different incidence angles 
and classified them into five different modes (Rossi 
et al. 2018, Durante et al. 2020). These results are 
also compared with small amplitude pitching 
NACA 0012 airfoil (Kurtulus, 2018; Kurtulus, 
2019). Rossi et al. (2018) and Durante et al. (2020) 
show that the classification of the wake structures 
behind NACA 0010 airfoil are similarly visualized 
for Reynolds number ranging from 100 to 3000 and 
at different incidence angles. 

In the present study, the focus is given to low 
angles of attack to understand the trailing edge 
separation phenomena at low Reynolds numbers of 
1000 to 4000. The novelty of the current paper is to 
show the effect of the Reynolds number at these 
low values where the data is not available in the 
literature to the best knowledge of the author on the 
instantaneous and mean flow fields and the 
aerodynamic coefficients for angles of attack up to 
10°. 


2. METHODOLOGY 


Unsteady, incompressible, laminar and two 
dimensional Navier-Stokes equations are used for 
the numerical simulations as is given by Eq. (1) and 
Eq. (2). A commercial software (ANSYS Fluent 
v14.0) which is a finite volume solver is used to 
obtain current results. SIMPLE-type fully implicit 
algorithm is used with pressure-velocity coupling 
and all the time and spatial solutions are selected to 
be second order accurate. 


V-V-0 (1) 


OV. es da - 
— -(V:V)V = --Vp + vV?V 
d í Q) 


The symmetric airfoils are pivoted from the 
quarter chord location and they are rotated with the 
corresponding angle of attack before the grid 
domain is formed. The field is composed of two 
domains, one close to the airfoil and the second one 
outside of the inner domain which extends 19c 
away of the airfoil towards both upstream and 
downstream (Fig. 1). 
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Table 1 Grid and time refinement studies for NACA 0012 airfoil at = 


Mesh Laminar/ Nodes Number of At [s] First cell Ca 
turbulent around elements spacing of BL 
analysis airfoil 


| 150 — 


(1994) 





With the preprocessor program, the airfoil is 
rotated for the given angle of attack in the inner 
region by keeping the outer domain and wake 
region fixed with the same structured grid. Three 
different meshes are compared for 10°, 20°, 35? and 
40? angles of attack in previous studies of the 
author and the results with medium and fine meshes 
are found to be very close to each other and also to 
the numerical simulations available in the literature 
(Kunz et al., 2000; Hoarau et al., 2003; Mittal et al., 
1994; Liu et al., 2012; Khalid et al., 2012; Suzuki et 
al., 2009). The boundary layer has a first spacing of 
0.002c for the medium mesh configuration selected 
whose number of elements are in the order of 2x10°. 
The number of nodes around the airfoil is 300 for 
the medium mesh configuration (Fig. 1). 

A time refinement study has also been 
performed (Kurtulus, 2015, 2016) and a time 
increment of At=0.005s is selected for the 
simulations. The average quantities are obtained in 
the interval of 50s <t <100s where quasi-steady 
solutions are obtained. Mittal et al. (1994) used for 
Reynolds numbers 1000 and 5000, a mesh with 
23600 elements and 23900 nodes around NACA 
0012 airfoil. Their unsteady numerical simulation 
method is based on the  Deforming-Spatial- 
Domain/Stabilized-Space-Time finite element 
formulation. The mean values of Mittal et al. (1994) 
in comparison with the current study grid 
refinement results are also tabulated in Table 1. 
Their results are found to be very similar with the 
results obtained in the current study. In the current 
study, fine mesh is also used for the highest 
Reynolds number and highest angle of attack 
investigated (Re=4000 and a=10°). The mean 
aerodynamic force coefficients are found to be very 
close to the medium mesh results as shown also in 
Table 1. 

The comparison of the mean aerodynamic force 
coefficients with the data available in literature 
(Mittal et al., 1994; Liu et al., 2012; Khalid et al. 
2012; Sunada et al., 1997) are given in Fig. 2. Mean 
lift coefficient for NACA 0012 at Re=1000 give 
very similar results with data available in literature 
(Liu et al., 2012; Khalid et al. 2012). The results at 


Re=4000 are also compared with the experimental 
data of Sunada et al. (1997) which are also available 
in Cleaver et al. (2012). Sunada et al. (1997), using 
a load cell, examined experimentally in a water tank 
the aerodynamic forces of fifteen different airfoil 
shapes one of them being NACA 0012 airfoil at 
Re=4000. Sunada et al. (1997) also observed that 
the roughness of the airfoils tested during the 
experiments does not affect the aerodynamic 
performance of the airfoils and this could be 
explained by the fact that this Reynolds number is 
far below the transition Reynolds number as is 
emphasized by Schmitz (1967). 


—O— Re=1000 

—@— Re=2000 

—®— Re-4000 

—vV-— Re=1000, Mittal et al. (1994) 

—*— Re=5000, Mittal et al. (1994) 
O Re=1000, Liu et al. (2012) 
<  Re-1000, Khalid et al. (2012) 
8G Re=4000, Sunada et al. (1997) 





0 2 4 6 8 10 
a [] 
Fig. 2 Mean lift coefficients in comparison 
with the literature data. 


The transition is defined as the passage from 
laminar and separated (subcritical) flow on the 
upper surface boundary layer of the airfoil to the 
turbulent and attached flow (supercritical). Schmitz 
(1967) has performed wind tunnel test with force 
measurements on the Reynolds number ranges of 
21000 to 168000 using 5 different airfoils. He 
denoted that for Re«150000, the airfoil has a critical 
range where the boundary layer become turbulent 
and the previously separated laminar boundary layer 
attaches. This situation results a sudden increase in 
lift coefficient and a decrease in the drag coefficient 
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and a high C/Ca ratio compared to Reynolds 
number regimes lower than transition. Reynolds 
number. 

The transition Reynolds number for NACA N60 
airfoil at a=6° was found to be at Re=63000 where 
Ci/Ca ratio in the supercritical condition becomes 
three times that of the subcritical case. From the 
stagnation point over the upper surface the airfoil, a 
negative pressure region is formed which is the 
main lift contribution region. From the minimum 
pressure region until the trailing edge, the pressure 
increases on the upper surface. The deceleration of 
the flow and the pressure increase can cause 
separation of the flow on the upper surface. On the 
lower surface of the airfoil, from the stagnation 
point to the trailing edge, and on the upper surface 
of the airfoil from the stagnation point to the 
pressure minimum, there is a tendency to maintain 
the laminar boundary flow due to the pressure 
increase and accelerated flow (Schmitz, 1967). For 
higher Reynolds numbers, above transition 
Reynolds number, from the minimum pressure 
location until the trailing edge, because of the 
pressure increase, there is also a tendency to 
transition to turbulence (Schmitz, 1967). 

A previous numerical study for lower thickness 
to chord ratios have been performed by Kunz et al. 
(2000) with laminar, two-dimensional 
incompressible Navier-Stokes solver for Reynolds 
numbers from 1000 to 12000. AII their numerical 
calculations use a C-grid topology with either 256 
by 64 cells or 512 by 128 cells. They also denoted 
that if there is no separation in the flow and the 
flowfield will be entirely laminar at these Reynolds 
numbers. Even from slight to moderate separation 
for Re«10000 will result a laminar reattachment 
(Kunz et al., 2000). 

The laminar results of Re=4000 are also 
compared with k-œ SST turbulence model with low 
Reynolds number correction using a low turbulence 
intensity of 0.1% for angles of attack of 0° to 10? 
and the aerodynamic force coefficients are found to 
be very close to the laminar simulations for this 
Reynolds number (Fig. 3). The laminar and 
turbulent analysis are also compared using fine 
mesh for the maximum angle of attack studied in 
the current study (a=10°) as tabulated in Table 1 
and results are also found to be close to each other. 

The mean pressure coefficient and the mean 
skin friction coefficient along the x direction are 
also shown in Fig. 4 for Re=4000 and maximum 


angle of attack studied in the current study (a=10°). 
The medium mesh and fine mesh solutions are 
Observed to be very close to each other. The 
distributions are also compared for laminar and 
turbulent solutions. The difference between the 
mean Cp and mean Cf coefficients all along the 
airfoil are also presented in Fig. 4c and Fig. 4d, 
respectively. The main differences are observed to 
be at the leading edge and then at the trailing edge 


of the airfoil. 
14 


0.8 








Fig. 3 
coefficient (b) at Re=4000 comparing laminar 
and k-o SST turbulence model (Tuz0.1 9c) 
solutions. 


Mean lift coefficient (a) and mean drag 
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Mean pressure coefficient (a), mean skin friction coefficient (b) and the discrepancies between 


two models ^C» (c), AC; (d) at Re=4000 and a=10° comparing medium mesh, fine mesh results with 
laminar solution and k-@ SST turbulence model (Tu=0.1%) solution. 


3. RESULTS 


The aerodynamic forces are highly unsteady 
depending on the separation type observed at the 
suction side of the airfoil. A detailed analysis is 
needed to understand the vortex shedding 
phenomena and oscillations of the aerodynamic 
forces at low Reynolds numbers. The amplitude of 
oscillation for the aerodynamic coefficients are 
found to increase rapidly with the angle of attack in 
these Reynolds number intervals. 

The minimum and maximum values of the 
unsteady lift and drag coefficients are presented in 
Fig. 5 with error bars in addition to the mean 
values. The oscillations become clear for Re=4000 
at critical angle of attack of Gc=4° and for Re=2000 
at acr =6°. The wake shedding is found to be Ocr=8° 
for Re=1000 for NACA 0012 airfoil. The vortex 
shedding phenomena is highly dependent to angle 
of attack, Re number and airfoil thickness (Kunz et 
al., 2000, Kurtulus, 2015, 2016; Rossi et al., 2018; 
Durante et al., 2020). The present study shows the 
effect of angle of attack and Re number on the onset 
of these oscillatory conditions. It is found that the 


lift curve slope at a=0° decreases with Reynolds 
number increase from 1000 to 4000 (Fig. 5a). The 
lift coefficient decreases as the Reynolds number 
increases from 1000 to 4000 for a<4°. The lift 
coefficient is found to be almost same for a=7° for 
these Reynolds numbers. For a>7°, the lift 
coefficients are found to increase with the increase 
of the Reynolds number. 

For the drag coefficient, there is a decrease of 
the mean values as the Reynolds number increases 
from 1000 to 4000 as is shown in Fig. 5b until 9? 
angle of attack. At a=9° angle of attack, Re=4000 
results higher drag coefficient with respect to 
Re=2000. At a=10° angle of attack, the increase of 
the drag coefficient is sharper and results higher 
values for Re=4000 compared to both Re=2000 and 
Re=1000 cases. As a result, a sudden jump of the 
lift to drag ratio is observed at a=10°. It should be 
noted that the turbulent solutions above 10? is found 
to change the aerodynamic force coefficients so this 
sudden increase in lift and drag coefficients are 
found to be as a precursor for the turbulent effects 
in 2D flows. It is also noted that at this angle of 
attack the amplitude of oscillations for both lift and 
drag coefficients are relatively big compared to 
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smaller angles of attack. This phenomena need 
further investigations therefore the current study is 
limited for the results where the laminar solutions 
are very similar with the turbulent results (see Fig. 3 
and Table 1). For all Reynolds number cases 
investigated, the lift to drag ratio of Re=4000 is 
found to be the biggest in the angle of attack 
interval of a.e[0? 10°]. Fig. 5c shows the Strouhal 
number based on the frequency (f) corresponding to 
the maximum amplitude of the amplitude spectrum 
of the unsteady lift coefficient obtained from the 
Fast Fourier Transformation as a function of the 
angle of attack and Reynolds number. The results 
are presented for both below the critical angle of 
attack (hollow symbols and dashed lines in Fig. 5c) 
and above the critical angle of attack (filled symbols 


0,25 4 
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and solid lines in Fig. 5c). The critical angle of 
attack is defined as the angle where the onset of 
oscillation is observed and a very major 
fundamental frequency can be detected from the 
amplitude spectrum of the aerodynamic force 
coefficients who shows oscillatory behaviors. A 
more detailed discussion for the critical angle of 
attack is given in the following section. Below the 
critical angle of attack, there is very low scale 
frequency amplitudes which does not show a 
dominant frequency. Some of these low scale 
frequencies below the critical angle of attack could 
be due to the internal corrections of the commercial 
solver aimed at reaching the convergence of the 
computations. 
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The mean lift coefficient (a), the mean drag coefficient (b) and Strouhal number (c) for NACA 


0012 airfoil at low Reynolds numbers. 


In the current study, the maximum frequency 
amplitude is also selected for the angles of attack 
below the critical angle to obtain Strouhal number 
where the wakes are not even shed and the 
aerodynamic force coefficients are not in a clear 
oscillatory behavior (see Fig. 15a for amplitude 
spectrum of lift coefficient at a=3°). For these low 
angles of attack a persistent laminar separation 
bubble (LSB) is observed. 

The mean lift and drag coefficient variations 
with respect to Reynolds number are also presented 
in Fig. 6. It is clearer from Fig. 6 that as Reynolds 
number increases, the lift coefficient decreases for 
ax4?. For 5°<a<6°, the minimum mean lift 
coefficient 1s obtained at Re=2000, so the mean lift 
coefficient decrease and then increase in the Re 
number range of 1000 to 4000. The mean lift 
coefficient is found to be almost constant for a=7° 
for Reynolds numbers investigated. For a28°, the 
mean lift coefficient is found to increase with the 
increase of Re number. A very sharp increase in 
mean lift coefficient is observed for a=10°. The 
mean drag coefficient is found to decrease with the 
increase of the Reynolds number for the angles of 
attack less than equal to 8°. For a=9° and a=10°, 
the minimum mean drag coefficient is obtained at 
Re=2000. As the Reynolds number increases both 
the lift and drag coefficient amplitudes are observed 
to grow above the critical angles of attack where 
oscillations exist. 





Q 1000 2000 3000 4000 5000 
Re a) 





0 1000 2000 3000 4000 5000 


Re b) 


Fig.6 | Mean Ci (a) and mean Ca (b) of NACA 
0012 airfoil at different incidence angles for 
Re=1000, 2000 and 4000. 
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3.1 LSB Formation Angle and Critical 
Angle for the Onset of Oscillations 


Oscillation in the aerodynamic coefficients are 
mainly dependent to the pressure and skin friction 
distributions on the airfoil which also changes 
instantaneously by the formation of the laminar 
separation bubble (LSB). At very low angles of 
attack, LSB is found to be close to the trailing edge 
and an increase in the angle of attack causes the 
separation point to move forward towards the 
leading edge (Fig. 7-Fig. 9). As also denoted by 
Tani (1964) a further increase of the angle of attack 
moves the separation point so far forward that the 
flow no longer reattaches to the surface within a 
short distance therefore the bubble has broken down 
or burst. Fig. 7 to Fig. 9 show the instantaneous 
pressure coefficient contours and streamlines at 
different angles of attack from Re-1000 to 
Re=4000 at t=100s. Below the critical angle of 
attack, it was also shown that the instantaneous and 
mean velocity fields reveals same vortex patterns 
(Kurtulus, 2015) as is also shown in Fig. 7. 

It is clear from Fig. 7-Fig. 9 that, for a=3°, LSB 
at the trailing edge is generated and grow as the Re 
number increases from 1000 to 4000. In addition, 
its separation point moves towards the leading edge 
as the Reynolds number increases. Both the mean 
lift coefficient and the mean drag coefficient are 


t*100 s (Re1000, 3 deg 


t5 1M s (Re 1000, 4 de 


Fig. 7 





found to decrease with the increase of the Reynolds 
number from 1000 to 4000 at a=3° as is also shown 
in Fig. 6. Although the loss of the mean lift 
coefficient is not very effective, the drag coefficient 
is found to sharply decrease even in the presence of 
a growing LSB at the trailing edge. Fig. 10 shows 
the instantaneous and mean pressure and skin 
friction coefficient distributions around the NACA 
0012 airfoil at a=2° for the Reynolds number range 
of 1000 to 4000. The instantaneous and mean 
distributions are found to be very close to each 
other. Skin friction coefficient all along the airfoil 
surface is found to decrease by the increase of the 
Reynolds number from 1000 to 4000 (Fig. 10b). No 
laminar separation bubble is observed at Re=1000 
for a@=2°. However when the Reynolds number 
increases to 2000, LSB is formed at the trailing 
edge of the airfoil. Trailing edge pressure 
distribution is found to decrease in magnitude as the 
Reynolds number increases (Fig. 10a). However, 
the suction pressure at the leading edge is increased 
in magnitude as the Reynolds number increases. 
The lift coefficient slightly decreases with the 
Reynolds number increase due to the decrease of 
the discrepancy between the upper and lower 
surface pressure distributions however, the skin 
friction coefficient reveals an effective drop of the 
drag coefficient with the Reynolds number increase. 


t=100.0 s (Re=1000, 6 deg) 


Instantaneous (t=100s, 1** row) and mean (2™ row) pressure coefficient contours and 


streamlines for Re=1000 at 3°<a<6°. 
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Instantaneous (t=100s) pressure coefficient contours and streamlines for Re=2000 at 0°<a<3°. 
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It is also very well-known that at these low 
Reynolds numbers separation starts from the 
trailing edge and as the angle of attack increases the 
separation point moves towards the leading edge 
creating a bigger separation bubble (Kunz et al. 
2000; Kurtulus, 2015, Eljack et al. 2021). 

LSB is first distinguished at an incidence of 4° 
for Re=1000 and found to persists up to higher 
angles of attack (Fig. 7). This angle is referred as 
LSB formation angle in the current study (ass). At 
higher angles of attack, for each Reynolds number 
and airfoil thicknesses, the LSB separation point 
moves forward towards the leading edge then turns 
to a leading edge separation which covers the whole 
airfoil upper surface by the formation also of the 
secondary, tertiary vortices. This onset of 
oscillatory behavior in unsteady force coefficients 
and alternating vortex shedding is referred as 
critical angle (Qcr) in the current study and the 





Instantaneous (t=100s) pressure coefficient contours and streamlines for Re=4000 at 0°<a<3°. 
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Fig. 10 Instantaneous (t=100s) and mean a) Cp and b) Ct distributions around NACA 0012 airfoil at 
a=2° for different Reynolds numbers 


shedding frequency at this angle of attack is denoted 
to be the fundamental frequency, fcr, at the critical 
angle of attack. The LSB is observed above an LSB 
formation angle and it is found to be stable and 
persistent characteristic for angles of attack lower 
than the critical angle. Just above the critical angle, 
the vortices are observed to shed alternatingly and 
the aerodynamic forces above this angle of attack 
have almost sinusoidal patterns for the alternating 
vortex shedding mode conditions (see Fig. 12b and 
Fig. 14). 

In a previous study of Kurtulus (2016), the 
vortex shedding patterns were classified into five 
different vortex shedding modes according to the 
vortex pattern, amplitude spectrum of Ci, and the 
longitudinal and lateral spacing of the vortices 
behind NACA 0002 and NACA 0012 airfoils at 
Re=1000 considering different angles of attack. 
Mode 1 which was denoted as ‘‘continuous vortex 
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sheet mode’’ is obtained at a<8° for NACA 0002 
and a<7° for NACA 0012. Mode 2 nominated as 
"alternating vortex shedding mode’’ was observed 
at 9°<a<11°, 14°<a<20° for NACA 0002 and 
8°<a<22° NACA 0012. It is thus noted that the 
critical angle also represents the angle when this 
vortex shedding mode changes from Mode 1 to 
Mode 2. Huang et al. (2001) observed an attached 
flow regime in their experimental study where they 
noticed an attached flow at the suction surface of 
NACA 0012 airfoil without a separation or vortical 
structure at a=2° and Re=1200. 

Table 2 shows in summary the LSB formation 
angle, critical angle and fundamental frequency for 
Reynolds numbers between 1000 and 4000 for the 
case investigated in the current study. No vortex is 
detected on the NACA 0012 airfoil surface for a<3° 
for Re=1000 so the LSB formation angle is found to 
be 4° at Re=1000 for NACA 0012 airfoil (Fig. 7). 
However, LSB formation is detected at a=0° at 
Re=4000 for the same airfoil (Fig. 9). Since a=0°, 
is a Symmetric condition both for the upper surface 
and the lower surface of the airfoil, the same trailing 
edge vortex is observed. Fundamental frequency is 
found to increase with the Reynolds number from 
1.28 Hz for Re=1000 to 9.34 Hz for Re=4000. Both 
the LSB formation angle and critical angle 
decreases with the increase of the Reynolds number. 


Table 2 LSB formation angles and critical angles 
for NACA 0012 airfoil at low Reynolds numbers 


LSB formation 4? 2> 
angle (arse) 


Critical angle 4? 
(Ger) 
1.28 Hz 3.56 Hz 9.34 Hz 


Flow regions are compared with the regions 
defined by Huang et al. (1995) for NACA 0012 
airfoil at different Re numbers and angles of attack 
as shown in Fig. 11. Region A represents Attached 
Boundary Layer, Region B represents Separated 
Boundary Layer and Region C is denoted as the 
Boundary Layer completely separated region where 
the flow completely separates from the airfoil upper 
surface. From the numerical simulations given by 
Kunz and Kroo (2000), where they emphasized that 
they stopped the simulations when the steady 
convergence is not reached. From the simulations 
that they performed, it can be found that the critical 
angle is above 7.5° for Re=2000 and above 6° for 
Re=6000 for 8% thick NACA 0008 airfoil. They 
found critical angles almost same for NACA 0006 
airfoil. The critical angle is above 7° for Re=2000 
and above 5° for Re=6000 for NACA 0004. As they 
performed also simulations for 2% thick NACA 
0002 airfoil, they stopped their analysis at 5.5°, 
which shows that the critical angle is above 8° for 
Re=1000, above 5.5° for Re=2000 and above 4° for 
Re=6000 for NACA 0002 (Table 3). A similar 
result is also found by Kurtulus (2016) for NACA 
0002 airfoil at Re=1000 where the critical angle is 
detected in the range of ac €[8° 9°]. It can be seen 
that the critical angle decreases as the Reynolds 





number increases. It is clear that the unsteadiness is 
delayed at higher angles of attack by the decrease of 
the Reynolds number which is also a result found in 
the current study (see Table 3). The transition angle 
from the continuous mode to alternating mode 
shedding has also been studied by Gopalakrishnan 
Meena et al. (2018) for symmetric NACA profiles 
for different thicknesses with and without Gurney 
flap at Re=1000. They compared their results with 
Kurtulus (2016) and found very similar results for 
NACA 0012 at Re=1000. Their findings about 
critical angle of attack 1s also tabulated in Table 3. 
Gopalakrishnan Meena et al. (2018) also performed 
3D analysis and found that alternating shedding 
mode tend to remain two dimensional even with 
spanwise effects at these low angles of attack and 
Reynolds number. 

Thickness decrease has an increasing tendency 
in critical angle for Re=1000 (Gca=8° for NACA 
0012, 8°<aer<9° for NACA 0002 and 7?«a«c «8? 
for NACA 0000). However, as the Reynolds 
number increases, the critical angle is found to have 
a decreasing trend for 1296 thick NACA 0012 
airfoil (Table 3). For Re=2000, the airfoil thickness 
is found to affect the critical angle with an 
increasing trend from 12% airfoil (5°<ocr <6°) to 
6% airfoil thickness (Qcr>7.5°) then the critical 
angle has a decreasing trend towards the 4% 
(Ocr>7°) and 2% thick airfoils (Qter>5.5°). 


1.6 O Region A-Region B boundary 
(Huang et al. (1995)) 
1.4 [D Region B-Region C boundary 
(Huang et al. (1995)) 
12 o LSB formation angle 

O| >= Critical angle 
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*Region A: Attached Boundary Layer, Region B: 
Separated Boundary Layer, Region C: Boundary Layer 
Completely Separated defined by Huang et al. (1995) 


Fig. 11 Reynolds number versus angle of 
attack map of NACA 0012 airfoil at low 
Reynolds numbers. 
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Table 3 Critical angle at different Reynolds 
numbers and airfoil thicknesses in comparison 
with the literature data 
angle (E): Experimental 

(N): Numerical 


1000 0.18 Oy =6° Gopalakrishnan Meena 
et al. (2018) (N) 


1000 Kurtulus (2016) (N) 


Gopalakrishnan Meena 

et al. (2018) (N) 
1200 
2000 
3195 
4000 
2000 


6000 


1000 


2000 
6000 
2000 
6000 


pane 
et al. (2018) (N) 

[ 006 | a5" | Kunzet al, 2000) (N) | 

[ 006 | a6" | Kunz etal. (2000) (N) - 


Q8? Kunz et al. (2000) (N) 


0.00 
(Flat 
plate) 





1000 


3.2 Instantaneous Aerodynamic Force 
Coefficients and Flowfields 


Instantaneous lift and drag coefficients for 
different Reynolds numbers investigated are 
presented in Fig. 12 for a=0° and a=8°. At a=0°, 
apart the impulsive start interval, the unsteady lift 
and drag coefficients are found to be almost 
constant for Re=1000, 2000 and 4000 (Fig. 12a). As 
the angle of attack increases to a=8°, the sinusoidal 
oscillations of the aerodynamic force coefficients 
are clearly visible for Re=2000 and 4000 cases (Fig. 
12b). The oscillatory behavior of the aerodynamic 
coefficients is mainly due to the both 
counterclockwise (CCW) trailing edge vortex and 
clockwise (CW) laminar separation bubble 
generated on the upper surface and their alternating 
shedding towards the downstream of the airfoil. The 
critical angle is denoted where the unsteady 
aerodynamic force coefficients exhibit an obvious 
unsteady oscillatory motion with an alternating 
vortex at the wake of the airfoil in addition to an 
amplitude spectrum of lift coefficient showing a 
distinguishable, a high amplitude peak value (see 
Fig. 15). It is found that for Re=1000, the critical 
angle of attack is 8°. The oscillation magnitude of 
the lift coefficient is -0.2% smaller (minimum Ci), 
+0.4% bigger (maximum C1) than the mean value 
for Re=1000. These oscillations increase abruptly 
with Reynolds number to almost -10% and +11% of 


Tars Gopalakrishnan Meena 
et al. (2018) (N) 
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the mean value for Re=2000 and -32% and +39% 
for Re=4000. 

The attached laminar boundary layer, which is 
subject to an adverse pressure gradient, results a 
separation of the flow which creates a laminar 
separation bubble on the airfoil surface at these low 
Reynolds numbers. For Re=2000 and a=8° case, 
the instantaneous streamlines and pressure 
coefficient contours are shown close to the airfoil in 
Fig. 13 to reveal the alternating vortex behavior 
between the time interval [99.5s 100s]. The 
instantaneous force coefficients at this time interval 
are also presented with a close view in Fig. 14. At 
maximum lift, drag and lift to drag ratio (at t299.5s, 
Cimax=0.3935, Camax= 0.1211), LSB (CW) is 
observed to start to detach from the upper surface of 
the airfoil (Fig. 13). The reattachment point at the 
trailing edge is pushed downstream of the airfoil 
and create a saddle point just downstream of the 
airfoil’s trailing edge (Perry, 1975). High suction 
pressure region is observed on the overall upper 
surface of the airfoil. At t=99.5s, a new trailing 
edge vortex is also observed to be formed due to the 
shear layer which curves from the lower surface 
towards the upper surface and it grows at the upper 
surface of the airfoil just at the trailing edge by 
pushing downstream the previously generated 
vortex (see t= [99.525s 99.60s] interval pressure 
contours in Fig. 13). The minimum lift coefficient is 
observed at t= 99.625s as shown in Fig. 14, where 
the trailing edge vortex grow and the LSB detaches 
from the airfoil surface (Fig. 13). This is also a time 
which is very close to the minimum drag coefficient 
point (Fig. 14). 

At t= 99.625s, as shown with the streamlines 
in Fig. 13, LSB is found to shed at the wake of the 
airfoil with a trace of LSB left behind it on the 
airfoil upper surface. At t=99.65s, LSB once more 
reattach to the upper surface of the airfoil and both 
trailing edge vortex and LSB is observable (Fig. 
13). Both the lift and drag coefficients starts to 
increase from their minimum values until t=99.77s 
(Fig. 14). The reattachment point of the LSB shifts 
towards the trailing edge by resulting the trailing 
edge vortex (CCW) to be shed downstream as 
shown in Fig. 13 for t=99.75s. The LSB is now 
observed to settle on the upper surface of the airfoil 
with a clear reattachment at the trailing edge. With 
the periodic shedding of the LSB (CCW) and 
trailing edge vortices (CW) a sinusoidal pattern on 
the aerodynamic forces are observed and an 
alternating vortex pattern is seen behind the airfoil 
at the downstream. It should be noted that because 
of the periodicity, Cp contours and streamlines at 
t=99.7s and t=100 s are almost the same (Fig. 13- 
last column of 3" and 4" rows). 
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Fig. 12 Instantaneous lift and drag coefficients at different Reynolds numbers at a) a=0° and b) a=8°. 
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Fig. 13 Instantaneous streamlines and C, contours at Re=2000, a=8°. 
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Fig. 14 Instantaneous Ci, Ca, Ci/Cà values at Re=2000, a=8°. 


3.3 Amplitude spectrum of the lift 
coefficient and fundamental frequency 


The oscillations reveal a clear fundamental 
frequency in the amplitude spectra of the lift 
coefficients as can be seen from Fig. 15, specifically 
for angles greater than equal to critical angle of 
attack as was tabulated in Table 2. Above the 
critical angle, a clear fundamental frequency and its 
sub-harmonics are visible in the amplitude spectra 
of the lift coefficients (see Fig. 15 for a=4° at 
Re-4000 with maximum fundamental frequency 


fcr-9.34 Hz, see Fig. 15 for a=6° at Re=2000 with 
fer=3.56 Hz and see Fig. 15 for a=8° at Re=1000 
with fe=1.28 Hz). 

In the current study, more intension is also 
given to the frequency spectrum of the lift 
coefficient for the angles above LSB formation 
angle (œrsB) which is 4° for Re=1000, 2° for 
Re-2000 and 0? for Re-4000 (Table 2). The 
maximum amplitude frequency gives a very low 
amplitude frequency spectrum without a noticeable 
peak as shown in Fig. 15 (see for a=3°). 
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.15 Amplitude spectrum of Ci for NACA 0012 airfoil for 3°<a<8° for different Reynolds numbers. 
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4. CONCLUSION 


The flow behind NACA 0012 airfoil at low 
Reynolds numbers has been studied at low angles of 
attack for the Re=1000, 2000 and 4000. The 
transition from the steady condition to periodic 
force evolution has been revealed with a detailed 
flow field analysis. The critical angle from the 
steady convergence condition (Mode 1) to unsteady 
oscillatory behaviors (Mode 2) were previously 
obtained at a=9° for NACA 0002 and a=8° for 
NACA 0012 at Re=1000. The fundamental 
frequency of shedding for alternating vortices is 
evaluated for different Reynolds number at the 
critical angle of attack. Another solution obtained 
from the current study is that although the mean lift 
coefficient is not affected too much from the 
incremental change of the Reynolds number from 
1000 to 4000 for NACA 0012 airfoil, the mean drag 
coefficient is found to decrease as the Reynolds 
number increases for the angles of attack less than 
equal to 8° even in the presence of a growing LSB 
at the trailing edge of the airfoil. 
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